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ABSTRACT 

We study, by means of MHD simulations, the onset and evolution of fast reconnection via the 
“ideal” tearing mode within a collapsing current sheet at high Lundquist numbers {S ^ 10^). We 
first confirm that as the collapse proceeds, fast reconnection is triggered well before a Sweet-Parker 
type configuration can form: during the linear stage plasmoids rapidly grow in a few Alfven times 
when the predicted “ideal” tearing threshold is approached from above; after the linear phase 

of the initial instability, X-points collapse and reform nonlinearly. We show that these give rise to 
a hierarchy of tearing events repeating faster and faster on current sheets at ever smaller scales, 
corresponding to the triggering of “ideal” tearing at the renormalized Lundquist number. In resistive 
MHD this process should end with the formation of sub-critical (5 < 10^) Sweet Parker sheets at 
microscopic scales. We present a simple model describing the nonlinear recursive evolution which 
explains the timescale of the disruption of the initial sheet. 


1. INTRODUCTION 

Magnetic reconnection is thought to provide the path¬ 
way for energy release in solar flares and other phenom¬ 
ena where energy is accumulated in the magnetic fields 
and currents of ionized high temperature plasmas. What 
has remained difficult to understand has been the trig¬ 
gering and speed of the process itself, which, given the 
extremely large Lundquist (S) and Reynolds {R) num¬ 
bers, implies that currents must collapse to extremely 
small scales before anything close to realistic timescales 
are approached. 

Indeed, an active region in the solar corona has a typ¬ 
ical spatial scale L of about L ~ 10®cm, a magnetic 
field B ~ 50G, density p ~ 10®cm“^ and a tempera¬ 
ture T ~ 10®K, hence a macroscopic Lundquist number 
S ~ 10^^. A flare emits a total energy of about lO^^ergs 
on a typical timescale of a few minutes, which cannot be 
explained by simple magnetic held diffusion or by a re¬ 
connecting instability occurring on the macroscopic scale 
itself. It has been suggested then that the reconnection 
trigger might only be provided by kinetic effects, beyond 
the (resistive) magnetohydrodynamic description of the 
plasma, or, alternatively, by the fast plasmoid instability 
of Sweet-Parker current sheets (|Cassak &: Drake! (|2013lf 
and references therein). 

It has been predicted recently, using linear stability 
analyses, that the tearing mode instability should grow 
on ideal timescales (“ideal” tearing, or IT) once the 
inverse aspect ratio ajL of a current sheet reaches a 
scale ajL ^ preventing the formation of the 
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paradigmatic Sweet-Parker current sheet (SP) for which 
Ogp/L ^ about 150 times smaller than 

for S ~ 10^^ and even smaller for greater S. Even 
more interestingly, the (linear) asymptotic tearing insta¬ 
bility in thin current sheets at arbitrary aspect ratios 
predicts that the maximum growth rate normalized to 
the macroscopic Alfven time increases nonlinearly with 
the aspect ratio, 'fm'Tp, ~ {L/a)^^^ (see also eq. ([T]) be¬ 
low): this suggests that collapse to the critical threshold 
thickness may provide the trigger for fast reconnection, 
where “fast” he re is meant for dynamics occurring on the 


ideal timescale (IPucci & Velli 

12014 iTenerani et al.ll2015t 

ILandi et al.ll2015HDel Sarto 1 

2nTa. 


At the intermediate Lundquist numbers usual to simu¬ 
lations, say, S ~ 10^ — 10®, the presence of plasma flows 
into (inflows), and along (outflows) the current sheet af¬ 
fects the scaling of the critical threshold for IT by induc¬ 
ing the formation of more elongated current sheets, that 
is, of layers havin g inverse asp ect ratios ajL ^ 5“““, 
with Uc > 1/3 ijVelli I l2015aD . Therefore, with the 
Lundquist number not sufficiently large, the distinction 
between the IT scaling and the SP scaling might seem 
academic. Nevertheless, the IT current sheet instability, 
as typical for the tearing mode, produces a quasi-singular 
inner layer Si, whose inverse aspect ratio follows the SP 
scaling Si/L ~ Such inner layer is however non¬ 

stationary, with inflows and outflows increasing exponen¬ 
tially in time—their ratio incidentally does not follow the 
SP scaling, but rather Uinjuout ~ —together with 

the amplitude of the reconnecting field: seen from this 
inner layer, one may interpret the dyn amics in terms of 
the “embedded reconnection” scenario (iCassak fc Drake! 
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I2009D . leading to greater energy storage and more effi¬ 
cient dissipation at the extremely large S values of the 
solar corona. 

Here we study the onset and evolution of the tear¬ 
ing instability within a single collapsing current sheet by 
means of resistive MHD simulations at S' = 10®. In the 
following paragraphs we first show that the transition to 
a fast tearing mode instability takes place during the col¬ 
lapse when the predicted threshold inverse aspect ratio 
a/L S“^/® is reached. Secondly, we show that the sec¬ 
ondary current sheets formed nonlinearly give rise to re¬ 
cursive tearing instabilities at increasingly smaller scales 
and faster super-Alfvenic timescales, which may also be 
well described by the flow-modified “ideal” tearing cri¬ 
terion and instability. The nonlinear evolution leads to 
a complete disruption in a timescale estimated at 0.05 
percent of the macroscopic Alfven time once initial tear¬ 
ing is triggered, which we model here in a corrected 
version of the fractal reconnectioii scenario first pro¬ 
posed by (|Shibata fc Tanuiiial 1200111 . Though the two- 
dimensional instability of current sheets and its subse¬ 
quent nonlinear evolution within resistive MHD has been 
studied before, it is shown here for the first time that the 
IT instability scenario, appropriately modified by the ef¬ 
fect of the reconnection velocity flows, provides a quan¬ 
titative description of the various stages in the evolution 
of a reconnecting current sheet, including hierarchical 
secondary island formation and disruption. Our results 
also provide a co herent framework through which pre¬ 
vious simulations (iLoureiro et al. 1 120051: iLanenta |[20?)M 
iDaughton et al. Il2009l : iBattachariee et al.ll2009l l may be 

more completely understood. 

2. TRIGGER OF FAST RECONNECTION: SCENARIO 

Let us consider a (local) current sheet of inverse as¬ 
pect ratio a/L, as the one given in eq. ([2]), with Alfven 
speed Va and magnetic diffusivity r], in a high S = 
Lv aIv ^ 1, constant density plasma. We assume, as 
in (jUzdenskv fc Loureirol[2r)T41 . that the current sheet is 
collapsing via some external process on a timescale Tc 
a{da/dt)~^, which we take of the order of the ideal 
timescale, Tc ^ = LJva- We extend the linear analy¬ 

sis of the incompressible tearing instability to this time- 
dependent case by defining an instantaneous growth rate 
7 (/c,t), k being the wave number along the sheet, where 
explicit time-dependence is introduced by the evolving 
aspect ratio L/a itself. 

In the static case, two regimes describe the unsta¬ 
ble spectrum of a current sheet: the so called small A' 
regime, for wave lengths close to th e instability thresh¬ 
old (or const. (jFurth et al.llI96,^ . for fca < I for our 
eq uilibrium), and th e large A' regime (or non const.- 
ijj (iCoDpi et al.llI976ll . for ka <C I). These two regimes 
have wave vectors k that lie to the right and to the 
left of the fastest growing m ode km., km. < k < 1/a 
and Q < k < krn. respectively (IBattachariee et al.l[20091 : 
ILoureiro et al. 11201,11 : iDel Sarto II20I5I 1. The wave vec¬ 
tor and t he growth rate of th e fastest growing mode are 
given by (|Pucci &: Vellill201^ : 

~ ^mTA ~ . ( 1 ) 

Similarly to eqs. o, the growth rates in the small and 
large A' regimes also increase with the aspect ratio; how¬ 


ever, their values tend to zero in the asymptotic limit 
S' —> oo, and the full 7 (fc) dispersion rel ation becomes 
rapid ly peaked around km (cfr. Fig. 1 of (|Pucci fc Vellil 
120141 1. Therefore, both the small and large A' regimes 
can be neglected in the framework of fast reconnection, 
since the development of the instability is controlled by 
the evolution of the fastest growing mode described by 
eqs. dD). Take the current sheet given by eq. ([2]), for 
example: the unstable spectrum lies within the range 
‘2'^/L < k < 1/a. From this condition and the first of 
eqs. CD it follows that the regime described by eqs. CD 
always exists once kmL > 27r, i.e. once a/L < 0.2 5' 

The second of eqs. CD shows that the growth rate be¬ 
comes ideal when approaching the critical width a/L ^ 


5“^/® (much smaller than 5“^/®) from above: the tran¬ 
sition to “ideal” tearing during collapse therefore occurs 
on the fastest growing mode. In conclusion, a collapse 
with Tc > Ta naturally drives a sudden switch from a 
quasi-stable state—growth rate depending on a negative 
power of 5—to an ideally unstable one. Otherwise, the 
disruption of the current sheet at widths thicker tha n 
critical (as considered in (|Uzdenskv fc Loureirol 1201^ 1 
would require an infinite time for 5 —> oo. 

We now consider the example of an exponential col¬ 
lapse. Exponentially thinning current sheets are ob- 
serve d in simulations of solar a nd stellar coronal heat¬ 
ing (iRappazzo fc Parkeil 1201311. in 3D MHD tur bu- 


lence ( Brachet et al.l l20I3t iGrauer fc Marlianil [20^ of 
interest to a broad range of astrophysical phenomena, 
as well as in th e nonliney stag e of the classic tearing 
instability itself (|Ali et al.ll20I4fl . 


3. NUMERICAL SET UP 

We consider a plasma with homogeneous density po 
and pressure pq. The background magnetic field Bq de¬ 
scribes an exponentially shrinking current sheet and is 
given by 

Bo = Bq tanh [y/a{t)] x + Bq sech [y/a{t)] z, (2) 

where the half width a{t) is prescribed and parametrized 
in time by 

a{t) = ao exp"*/”” -|-aoo(l — exp"*/””). (3) 


We employ a 2.5D compressible MHD code periodic in x 
and with non-reflecting bou ndary conditions along the 
inhomogeneous y direction (iLandi et al.ll2005ll . We as¬ 
sume an adiabatic closure, with index F = 5/3, a scalar 
resistivity, and a Newtonian viscous stress tensor with 
Prandtl number P = I. With this choice , viscosity does 
not a ffect the scalings of IT significantly (jTenerani et al.l 
I20I5I1 . The collapse of the background magnetic field is 
obtained by adding a source term P of the form 





dBo 


( 4 ) 


in Faraday’s equation. Magnetic and velocity fields are 
normalized to Bq and to = Bq/^/Attpq, respectively, 
lengths to the macroscopic length L, time to Ta, and 
density and pressure to pq and Pg/dTr, respectively. It 
is useful to introduce also the normalized flux function 
-tjj, such that B = V X i/z. We set po = 0.8 and 5 = 
10®. Instability is seeded with a random noise of small 
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Table 1 

Simulation parameters: Lundquist number, collapsing time, 
asymptotic half width, onset time, nonlinear time, and half width 
reached at the end of the linear stage. 


run # S 

rcpf. 


r,lT^ 

'^nl /"^A 


run 1 

1 


3 

16 

0.01 

run 2 

4 

5-1/3 

11.5 

22 

0.0105 

run 3 

10 

5-1/3 

16 

36 

0.013 

run 4 

1 

5-1/2 

3 

4.5 

0.0024 


amplitude. Because of the wide range of scales to be 
resolved in the y direction, we employ an inhomogeneous 
grid with increasing resolution at the neutral line. The 
simulation box has normalized dimensions x Ly = 
27r X 0.96 (runs 1-3) and L^x Ly = 2'kx 0.46 (run 4), with 
2048 X 1024 mesh points. Resolution at the neutral line is 
Ay = 0.0001, which allows to resolve the diffusion region 
of th e IT, that scales as 5i/L ~ (jPucci fc Vellil 

l2014fl . Since for a ~ L the instability is extremely slow, 
we start from oq = 0.1 L. This is a good compromise for 
limiting computational time while retaining a sufficient 
dynamical range for the collapse. 

4. RESULTS 

In Table [T] we list the background and main dynamical 
parameters for each run, and in Fig. [T]we show for ref¬ 
erence the temporal evolution of some unstable Fourier 
modes of the flux function ipkiv)-, in y = 0, from run 2. 
The growth rate and wave number of the fastest growing 
mode of the primary tearing at the critical IT threshold 
will be labelled ■ji and ki, respectively. In particular, in 
our simulations we find kjL = 10 and TiT. = 0.46, in 
agreement with theory (|Tenerani et al.l[2015h . 

4.1. Linear stage 

Tearing onset takes place at time r* such that 
7 m(T*)Tc = 1, where 7m = 0.46 After 

onset, the most unstable modes grow according to the 
WKB uolution ipm(t) ~ exp[/ 'y,n(t')dt'], represented by 
the dashed line in Fig. [T] In runs 1-3 we consider differ¬ 
ent cases of collapse in which Coo/T = In all these 

runs the linear stage is ultimately dominated by modes 



time (xj 

Figure 1. Upper panel: a/L vs. time. Lower panel: temporal 
evolution of the amplitude of some unstable modes 'ipkiv) the 
neutral line y = 0 (run 2). 


close to the ideally unstable one kiL = 10. Run 4 forces a 
collapse of the initial sheet down to aoo/L = but 

proves that inverse aspect ratios ajL cannot 

be formed: the tearing mode indeed ignites right after 
t = 2.4 Ta, at which a/L = and a large number 

of islands rapidly pop up. Even at the end of the linear 
stage (see Section the collapsing sheet thickness is 
twi ce the corresponding SP thickness at the given S and 
P (|Tenerani et al.l[20T^ . 

4.2. Early nonlinear stage 

Nonlinearity becomes important at r„;, when the 
width w of magnetic islands is of the order of the 
width of the inner diffusion layer 5. The fastest grow¬ 
ing mode has 5m which scales with the aspect ratio 
as 5m/a ~ [avp,/ri)~^/^ (iLoureiro et al. I [201311 . hence 
5m/a ~ S~^^^{L/aY/^. By using the definition w/a 
2 y/ipm/a (jBLskampI 1200011 . we estimate the amplitude 

~ 0.25 5'“^/^ [a(T„j)/L]^^^ ~ 2.5 x 10“^, in 
agreement with simulations. No nlinear effects lead t o 
the competition of two processes (jMalara et al. |[T9M^ : 
on the one hand, islands start to merge through an in¬ 
verse cascade process from the fastest growing modes 
to smaller wave numbers, see, e.g., the inverse cascade 
from kiL = 10 towards kL = 5 and below starting at 
Tni = 22 Ta (Fig. [T|); on the other hand, X-point collapse, 
and subsequent secondary current sheet formation, takes 
place during the further nonlinear growth of the islands, 
as is typical for strongly unstable modes far from the 
small A' regime (iJemella et al.ll2003f). in lieu of the slow, 
algebraic Rutherford growth (lRutherfordlll973[l . 

4.3. Nonlinear stage: recursive X-point collapse 

We analyze the fully nonlinear stage using run 1 
as a reference, shown in Fig. [21 X-point col¬ 
lapse proceeds leading to the f ormation of sec o ndary 
plasmoids (ma gnetic islands) ([Loureiro et al. I 120051 : 
lLaDenta1l2008|l and, next, a recursive process of super- 
Alfven ic secondary layer formation and disruption, takes 
place (iShibata fc Taiimnall200H iDaughton et al. I I2009I 
iBattachariee et al.ll2009H . Fig. [2] shows a temporal se¬ 
quence of recursive plasmoid formation, taking place 
at the center, near the flow stagnation point of the 
original instability. Previously, such recursive recon¬ 
nection has been modeled as a succession of unstable 


SP layers (ILoureiro et al. 1 120051: IDaughton et al. 1 

2009 

Cassak & Drake 

120091: Huang & Battacharieel 

2010 

Uzdenskv et al.l 

20101: ILoureiro et al.l 1201211. Our 

simu- 


lations show instead that it is driven by onset of the IT 
mode, triggered by the dynamical lengthening of sheets 
to the local critical threshold, in a way similar to that 
discussed in section 14.11 As the recursive X-point col¬ 
lapse occurs in steps, we label the half length of the n-th 
current sheet the inverse aspect ratio a„/L„, the nor¬ 
malized (half) width of the inner diffusion layer 
the local Lundquist number = (L„/L) S and Alfven 
time Ta,„ = {Ln/L)Tf,. Referring to Fig. [21 the recursive 
X-point collapse starts with the formation of Li from the 
primary diffusion region of thickness 5i, hence oi ~ 5i. 
Current sheet Li becomes unstable, less than 3 Ta af¬ 
ter the end of the linear st age (second panel) , consistent 
with what was found in (|Landi et al.l 120151) . Here we 
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Figure 2. Blow-up of the recursive “ideal” tearing instability in 
run 1: contour plot of the out-of-plane current density jz (color 
coded) and magnetic field lines (white). Notice that the whole 
domain is much larger, 0 < x/L < 27r and —0.48 < y/L < 0.48. 


show that IT triggering within Li induces the growth of 
two plasmoids and of a second current L 2 , with a 2 — 
(third panel), which is itself destroyed by multiple plas¬ 
moids in about 0.5 Ta (fourth panel). In Fig. [3] we plot 
By{y,XQ), in Xq/L = 3.2, at three different times. The 
shown profiles of the magnetic field closely resemble the 
tearing mode eigenfunction, and display a clear hierarchi¬ 
cal structure within each of the inner diffusion layers. To 
further analyze the tearing onset inside such secondary 
layers, in Fig. |4]we display Li/L (in green), ai/L (ma- 



y/L 


Figure 3. Profile of By{y, xq) in xq/L = 3.2 at t = 18 Ta (primary 
tearing, n = 0), t = 19Ta (n = 1), and t = 19.8Ta (n = 2). The 
inset is a blow-up of By{y,X()) at t = 19Ta and t = 19.8Ta. 


genta), and aijLi (in blue) as a function of time, from 
t = 15 Ta, up to t = 19 Ta. During the collapse process the 
width remains almost constant, ai ~ <5^ ~ (here 

Si ~ 0.0015), while the length Li stretches at about a 
rate Li/Li ~ Tc^i — 2ta. Fig.|3]shows that ai/Li 

_ 1 /Q 

crosses the critical threshold ' (red dotted line) at 
t ~ IOta. At that time Li ~ 0.1 L, thus Ai ~ 10®, 
P = 1, and linear theory predicts 7 iprA,i — 0.47, or 
7 inTA — 2.3 ^ t~iTa. In the light of the previous discus¬ 
sion, we would expect Li to disrupt more rapidly than in 
~ 3 Ta, though on the same order of magnitude timescale. 
On the other hand, ther e is a stabilizing effe ct from out¬ 
flows forming along Li (IBulanov et al.lll978f) : these may 
induce the formation of thinner sheets having inverse as¬ 
pect ratio Un/Ln ^ with ac > 1/3. Based on 

the empirical observation that flows stabilize SP up to 
Sc — 10^, one can extend the IT theory to find a criti¬ 
cal ac which includes the effect of inflow-outflows (jVelli I 
I2ni5ail : 


2 log /r -I- log S„ 
3 log S„ 


(5) 


In eq. ([5]), /r = r„/(/ 7 i) is the ratio of the plasmoid 
evacuation rate to the maximum growth rate, corrected 
by a thresho l d fac tor / ~ 0.5 — 0.1. The value fi = 
10 (iBiskamd 1200011 yields the observed ac = 1/2 for 
Sn = 10^, while as expected lim ac = 1/3. Note that, 

Sn—^OO 

as shown below, Sn is a decreasing function of n. The 
black dotted and dash-dotted lines in Fig. SI correspond 
to 5']"“° at two different /r, while the light blue dotted line 
corresponds to agp.i/Bi. Though our Lundquist numbers 
are not extremely large, Li still disrupts before reaching 
the SP width. The trend of the data plotted in Fig. S] 
goes in the direction of our scenario, that is, that the 
fast tearing instability is triggered within collapsing cur¬ 
rent sheets once the critical threshold is reached. Tak¬ 
ing into account the increased Alfven speed due to pile- 
up lust outside the inner diffusion region as considered 
by (|Cassak fc Drakell2009ll does not change the trend of 
our results. Interestingly, the di scussion of flux pile-up 
and embedded reconnection led (jCassak fc Drakel 1200911 
to independently find a similar scaling ~ 5“^/® for the 
inverse aspect ratio of unstable current sheets starting 
from an initial secondary Sweet-Parker sheet. Their anal¬ 
ysis starts from a different line of thought, in the search 
for the proper criterion to destabilize an embedded SP 
sheet, and proceeds via a linearization of the fields in the 
neighborhood of the neutral line. It then finds a scaling 
with S for the aspect ratio of the embedded layer scal¬ 
ing in a way similar to IT. On the other hand, the IT 
scaling derives from requiring an S'-independent growth 
rate from the complete eigenmode analysis of tearing in¬ 
stability theory, and finds that the entire sheet inverse 
aspect ratio scales as The results therefore have 

different origin, and for the IT case the presence of flows 
modifies the scaling exponent ac according to eq. : in¬ 
cidentally, the latter expression is in excellent agreement 
with the aspect rati os at which plasmoids are observed 
to be ejected in the (iCassak fc Drakell2009ll paper. 


5. DISCUSSION 
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time (xa) 

Figure 4. Length LxjL (green), width aijL (magenta), 

and inverse aspect ratio aijLi (blue) of the first secondary 
current sheet vs. time (run 1). The red dotted line represents the 
critical threshold and the light blue dotted line the 

viscous asp/L. Black dotted and dash-dotted lines represent the 
critical threshold with outflow effects, given by eq. ©■ 


We have shown that a collapsing current sheet disrupts 
due to the triggering of fast reconnection as predicted by 
the IT model. This process is shown to proceed recur¬ 
sively, giving rise to a hierarchy of fast tearing instabil¬ 
ities at ever smaller scales. Let us model the recursive 
collapse of X-points, in the spirit of the original fractal 
reconnection model, with the following differences: first, 
we don’t consider SP as the initial condition; second, we 
use the width a„, derived from the IT scalings, to find 
Ln- Neglecting for simplicity viscous and outflow effects, 
imposing that a„/Tn-i (i.e., a„ — Sn-i), and 

_1 /q 

that IT onset occurs at an/Ln Sn ' , we hnd: 

Ln = Sn = (6) 

For n —)► oo then Ta,„ —>■ and, as expected, 

S'n —?► 1 and ttn/Ln —>■ I. However, after a number n* of 
steps Sn^ — 10^, at which we expect to reach a stable SP. 
For S = 10^^ (solar corona) n* ~ 4, and the length 
suddenly drops to microscopic scales, Li/L ~ 6 x 10“^, 
L 2 /L ~ 2 X 10-®, Ls/L ~ 3 X 10"®, L 4 /L ~ 10-9. 
The recursive X-point collapse after the first trigger thus 
lasts a time interval Ttot — Ln=i Ti.n — 5 x IQ-^Ta, 
one order of magnitude less than the upper limit given 
by (jShibata fc Tanumall200l(l . Higher resolution simula¬ 
tions are nevertheless necessary to assess how the non¬ 
linear evolution saturates. Our scenario can be modi¬ 
fied to include kinetic scales which might b e reached dy- 
namically in typical astrophysical systems (jCassak et al.l 


l2005t iDaughton et al. 1 1200911 . In this case we expect a 
change of the power laws (|6]) as different scalings de¬ 
scribe t he physics at the critical threshold in the kinetic 
regime (jPel Sarto 11201 , 'tH . 

We thank D. Del Sarto for interesting discussions and 
useful comments. This work was supported by NASA 
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